Asymmetric solitons in saturable nonlinear medium with complex potential 
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We study the stability of soliton waves through saturable nonlinear inhomogeneous media with 
gain and loss, not necessarily in balance. We report on the existence of families of stable asymmetric 
spatial solitons in a saturable nonlinear medium in the presence of a gain/loss asymmetrical profile. 
A tunable parameter defines the deviation from the TT-symmetric case, and a thorough study on 
the existence and stability of soliton families is carried out. The properties of the nonlinear modes 
bifurcating from the eigenvalue of the underlying linear problem are investigated. The eigenvalue 
ranges in the power-eigenvalue diagrams for different gain/loss profile are inspected. We find that 
the saturable nonlinearity severely restricts these ranges as the eigenvalues tend quite fast to an 
asymptotic profile, as power increases. The examples of the dynamics of the asymmetric solitons 
obtained by numerical simulations exhibit a good agreement with the analytical results. 

PACS numbers: 42.65.Tg, 42.65.Sf 


I. INTRODUCTION 

The past decade has witnessed an intense activity 
in linear optical systems, described by non-Hermitian 
Hamiltonians, triggered by reports on a wide class of 
complex potentials that obeys the parity-temporal sym¬ 
metry, the well known TT symmetry [1]. These TT- 
symmetric Hamiltonians exhibit real spectra even though 
non-Hermitian, as it has been demonstrated in semi¬ 
nal works within a quantum context [2, 3]. However, 
these Hamiltonians are quite suitable to describe TT- 
symmetric optical structures, through the association of 
the local potential function with complex-valued refrac¬ 
tive index profiles. This association has been quite fruit¬ 
ful and has turned the field of optics and photonics to be 
a particularly adequate ground to test the amazing con¬ 
sequences of TT-symmetric systems. Actually, TT sym¬ 
metry was investigated theoretically and firstly demon¬ 
strated experimentally in the optical context [4-7]. The 
inclusion of gain/loss in nonlinear propagation is an even 
more interesting problem, particularly, in the case of 
the nonlinear Schrodinger equation. Based on the non- 
integrable character of the nonlinear Schrodinger equa¬ 
tion, Wadati has demonstrated that a system where TT 
symmetry is modeled by a refractive index having a sym¬ 
metric real part and an anti-symmetric imaginary one, 
admits soliton propagation [8]. This notion developed 
into a wide class of TT- symmetric as well as non-TT- 
symmetric potentials with real spectra. In the linear the¬ 
ory such modes were studied for exactly integrable mod¬ 
els in [9, 10], and their nonlinear extension were reported 
in [11]. Similarly, these classes of potentials have been 
investigated within the context of the Gross-Pitaewskii 
equation [12]. Recent reviews on the matter may be 
found in Refs. [13-15]. Various types of nonlinearities 
have been investigated: solitons in TT-symmetric local¬ 
ized potentials with second-harmonic generation were in¬ 
vestigated in [16] and was found that such system admits 
stable solitons above the TT-symmetric phase breaking 


threshold; gap solitons in one-dimensional TT -symmetric 
periodic potentials with saturable nonlinearity have also 
been reported in Refs. [17, 18]. Furthermore, a logarith¬ 
mically saturable nonlinearity investigation indicated the 
existence of stable solitons in TT - symmetric periodic 
potentials [19]. Quite recently an investigation on the 
fundamental and higher-order families of solitons in one 
dimension, in the presence of a TT-symmetric potential 
and a saturable nonlinearity, concluded that by increas¬ 
ing the saturation parameter one may partly hinder the 
instability of these solitons [20, 21]. 

The developments reported in TT-symmetric optical 
structures have also triggered further studies to include 
non-TT-symmetric complex potential functions in vari¬ 
ous contexts widening the potentiality of features such 
as a real spectra as well as a phase transition [22-24]. 
Several classes of non-TT-symmetric complex potentials 
of the Wadati type with real spectra have been suggested 
[25]. These potentials contain free parameters that rep¬ 
resent tuning. By tuning the refractive index below a 
critical value the spectrum undergoes a complex phase 
transition. By adding a real part with an adjustable pa¬ 
rameter one opens an additional tunable possibility with 
ranges that admit real spectra and also a phase transi¬ 
tion to the complex phase. Nonlinear Kerr systems have 
been studied in a non-TT-symmetric potential and sta¬ 
ble soliton solutions were found [26]. In Ref. [27] the 
existence of continuous families of nonlinear modes in a 
complex asymmetric potential are attributed to an un¬ 
derlying conserved quantity in the system’s dynamics. 
Based on this fact, the existence and stability of solitons 
in quadratic nonlinear unidimensional media under the 
influence of a one-parameter family (an asymmetric class 
of Wadati potential) have also been studied [30]. Families 
of stable solitons bifurcating from the fundamental mode, 
as well as from the second harmonic, were reported. Ma¬ 
terials with high nonlinear coefficients such as semicon¬ 
ductor doped glasses and organic polymers, cannot be 
described by a Kerr-type nonlinearity [28], as in these 
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materials, the saturation of the nonlinear refractive in¬ 
dex occurs at moderately high pulse intensities, so that 
one must study the effect of higher order nonlinear it ies. 
Asymmetric solitons due to the symmetry breaking phase 
transition have been studied in a saturable medium [29]. 
However, to the best of our knowledge, a saturable non¬ 
linearity in non-TT-symmetric potentials have not been 
under investigation. Furthermore, theoretical and exper¬ 
imental results have indicated that a saturable nonlinear 
medium is suitable for the propagation and control of 
spatial solitons. As a notorious example photorefractive 
optical solitons, which can be induced at relatively low 
power formation, are quite important for applications in 
all-optical switching and signal processing [31]. There¬ 
fore, in this paper we propose an asymmetric potential 
to describe the existence and the stability properties of 
solitons in a saturable nonlinear medium. Our main ob¬ 
jective in this work is to include an asymmetry parameter 
to see whether one may attain a flexibility with respect 
to the balance between nonlinearity and asymmetry that 
stabilizes the soliton. In this way we hope to extend 
the class of potential functions, manipulating, satura¬ 
tion, gain and loss, to achieve new functionalities. This 
paper is organized as follows. In Section II the theoret¬ 
ical model is described. Section III briefly presents the 
linear properties of the potential function. The stability 
and dynamics of the localized solutions are considered in 
Section VI. 


II. THEORETICAL MODEL 

Let us begin by writing the propagation equation of an 
optical field ?/(£,£) through a diffractive nonlinear sat¬ 
urable medium: 

.du d 2 u r 9 . -i Id 2 ^ \ 

l d( + W + ^ 9 +a9 + U + a i + 7 | m |2 u = ° ^ 

u being the dimensionless amplitude of the field while 
£ and £ are, respectively, dimensionless transverse and 
propagation coordinates, scaled to the characteristic size 
of a localized modulation of the refractive index, which is 
fully described by a real function g = #(£). The real part 
being g 2 +ag , where a is a free parameter which plays the 
role of a relative strength of the real part of the poten¬ 
tial to the imaginary part, which is defined as: g% = ^|. 
This adjustable parameter <a, widens the scope of non- 
Hermitian Hamiltonians with real spectra and may drive 
the system to a phase transition to the complex plane 
[25]. Fig. 1 shows the variation of the real part g 2 + ag 
as a function of a illustrating its effect on the amplitude 
of the real part of the potential function. The last term in 
equation ( 1 ) represents a saturable third order nonlinear 
material with saturation parameter 7. The case 7 = 0 
represents the usual Kerr nonlinearity without satura¬ 
tion. The nonlinearity may be of a self-focusing (a > 0) 
or a self-defocusing type (a < 0). 


We look for localized solutions of (1) in the form, 

«(£,C)=™(0e i6C ,, (2) 

where b is the soliton propagation constant and w(£) is 
the solution of 

d 2 w r 9 7l \w\ 2 , . 

df 2+19 +ag + ig i -b\w + a ^- : -^ w, (3) 

constrained to zero boundary conditions, that is, w (£) —)> 
0 as |£| —>> oc. 


III. LINEAR PROPERTIES OF THE 
ASYMMETRIC LOCALIZED POTENTIAL 

Previous work [25] has shown that, similarly to the 
TT-symmetric case, the complex potential 

G(£, a) = g 2 + ag + ig% (4) 

when substituted into the linear version of equation (3), 

Lw = biw, (5) 

where b\ is the linear eigenvalue and the operator L is 
defined as 


L ~^+°, ( 6 ) 

depending on the value of ck, may exhibit all-real spec¬ 
trum. Furthermore it was proved that, for sufficiently 
mild conditions the value a = 0 always guarantees a real 
spectrum for ( 6 ). Let us now briefly describe the recently 
reported [30] linear properties of ( 6 ) when the potential 
G is generated by the following function: 


g = V sech 




( 7 ) 


where the constant V is the amplitude and the parame¬ 
ter (3 represents its degree of asymmetry. The function g 
is continuously differentiable, has a global maximum at 
£ = 0 and is asymmetric for /3 ^ 0 , as is illustrated in 
Fig. 1. The function g decays as when £ —>• — oo 
and when £ —>> oo. In Fig. 1 we show the 

spectrum of ( 6 ) using (7) for different values of a and 
a fixed f3 = — 0.5. A complex phase transition is found 
for a < a cr = —0.187, where a pair of complex eigen¬ 
values emerge due to the collision of two localized eigen¬ 
values at bi = 0.013 (See Fig. 1 ). Note that for a > 0 
there are no phase transitions, just the appearance of 
additional localized modes due to the higher amplitude 
of G. Here, in contrast with the symmetric system for 
which /3 = 0 , there are no degenerated discrete eigenval¬ 
ues as there is no symmetry to induce degeneracy. In¬ 
stead, as a changes continuously, a complex phase may 
appear only due to the collision of simple discrete modes 
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FIG. 2: Discrete spectrum of b\ as a function of f3. Note the 
appeareance of a continuous band at /3 = —1. The values 
chosen for the potential are V = V0.5 and a = 0. 


IV. NONLINEAR LOCALIZED MODES AND 
LINEAR STABILITY ANALYSIS 


FIG. 1: Upper left panel shows the eigenvalues of (6) with 
V = 0.5, ^ = —0.5 as a function of the parameter a. Up¬ 
per right panel shows the effect of a on the shape of the real 
part of the potential G. Lower panels show the eigenvalues 
below (a = —0.18) and above (a = —0.2) the phase transi¬ 
tion, where the real spectrum breaks into a complex phase. 
Note the appearance of a complex pair of eigenvalues due to 
the collision of two real valued eigenvalues at b = 0.013 where 
a cr — —0.187 (not shown). The spectrum contains an imagi¬ 
nary pair for a < a cr - 


as illustrated in Fig. 1 . Albeit shown here for the spe¬ 
cific case of f3 = —0.5, different values of (3 are found 
to exhibit the same behavior with respect to the break¬ 
ing of the real valued phase of the spectrum. The effect 
of different values of the asymmetry parameter /? on the 
eigenvalues b\ is shown in Fig. 2. In all results presented 
below, we have chosen the values V = dzy/Odb Note in 
(7) that this means that the absolute value of the am¬ 
plitude of g is fixed, so that only the decay rate of the 
£ > 0 region is modified. This affects the values of b\: 
as (3 decreases, the area of g becomes larger and, conse¬ 
quently, the number of localized modes increases. Clearly 
if /3 = —1, then 


lim G(£) = 0 , 

( 8 a) 

— 00 

lim G(0 =V 2 + aV. 

( 8 b) 

£—>-00 

This is an interesting situation as for (3 = 

— 1 the poten- 


tial G becomes very similar to a Heaviside step function 
with some localized imaginary part while the spectrum 
approaches a continuous band defined by b\ < V 2 + aV. 


Let us now, investigate the existence and stability of 
the localized solutions of (3). From now on we restrict 
our considerations mainly to solutions bifurcating from 
the linear limit, more specifically, as w 0 , b tends to a 
finite eigenmode of the linear problem i.e., b b\. The 
total energy P of a localized solution is defined by: 

POO 

P= / \w\ 2 <%. (9) 

J oo 

Let us now perturb the solution ?/(£, £) by writing: 

u&O = (M0+P+(0e- iA< +P-(0e iJ< ) e iK , (10) 

where, |p±(£)| 1, represents small perturbations. In¬ 

serting ( 10 ) into ( 1 ) and collecting linear terms in p±(£) 
one ends up with an eigenvalue problem given by: 



L s 


L — 


b + a 


M 2 (tM 2 +2) 

(1+7M 2 ) 2 


—a 


(i+ 7 M 2 ) 2 


—L 


u (i+ 7 M 2 ) 2 

7 _ M 2 (7M 2 +2) 

( 1 + 7 | k ;| 2 ) 2 


(lib) 

so that, whenever an eigenvalue A with Im(A) > 0 occurs 
the solution is linearly unstable. 


A. Effects of saturation parameter 7 

Investigating the effect of the saturation parameter 7 
we have found that all investigated branches with 7 > 0 
exhibit a vertical asymptote in the ( 6 , P) plane as shown 
in Fig. 3 and Fig. 4. 

This is expected since for solutions with very high in¬ 
tensities one has, 


lim 

\w(£)\ 2 ^oo 


HOI 2 

1 + 7KOI 2 


a 

7 


( 12 ) 
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FIG. 3: Total energy P as a function of b with V = y/oH, P = 
—0.5, a = 1 and different values of the saturation constant 7 . 
Dashed lines represents the asymptotes 61 +< 7 / 7 . All solutions 
are stable. 




FIG. 4: Total energy P as a function of b with V = a/T5, P = 
—0.5, cr = — 1 and various values of the saturation constant 
7 . All solutions are stable. 


Now, inserting (12) in (3) we obtain 


cPw 

~de 


+ 


g 2 + ag + ig^ - b+- 


w, 


which, of course the same as (5) 




Lw = (6-)re 

7 


(13) 


(14) 


and has a localized mode at b = b\ + where b\ is a 
discrete value of (5). In other words, as the amplitude 
of the soliton increases its profile becomes similar to the 
shape of the linear eigenmode it bifurcated from, in the 
limit re —>• 0. 

In fact, extensive numerical calculations of branches 
show that the fundamental branches exist in the re¬ 
gion b G (61,61 +7 -1 ) for the a = 1 case and 6 G 
(61 — 7 _1 ,6 i) when a = —1. Interestingly enough as, 
using the fact that the branch must be localized within 
the region between 61 and &i + -, one realizes that the 
a = l(cr = — 1) branch bifurcates to the right (left). 

In what concerns stability, we have investigated two 
main potential functions, i.e., G(£) and G*(£) which cor¬ 
respond to the cases V = \/043 and V = — \/0~5 respec¬ 
tively, if one also substitutes a —a. It should be noted 


FIG. 5: Total energy P as a function of b with V = — v / T5, 
P — — 0.5, cr — 1 and for various values of saturation constant 
7 . Right panel shows an expanded view of the region where 
the transition between stable and unstable solutions occurs. 
Thick lines represent linearly stable solutions and lines repre¬ 
sent unstable solutions. 


V = WO.o, a = 0, 3 = -0.5, a = 1 



FIG. 6 : Maximum value of imaginary part of (11) as a func¬ 
tion of the total energy of localized modes P for several values 
of 7 . Note that all branches are unstable but max{ImA} —> A 7 
as P —>> 00 . The higher the value of 7 , the lower the value of 
A 7 . Parameters are V — P — —0.5, a = 1. 


V = -V0.5, a = 0,6= -0.5, a = -l 



FIG. 7: Total energy P as a function of b with V = — v / R5, 
P — —0.5, <7 — —1 and different values of saturation constant 
7 . All solutions are stable. 


here that, if re is a solution of (3) with G then clearly re* 
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is a solution of (3) with G*, and both solutions are ob¬ 
tained for the same value of P. Despite no difference in 
the (6, P)— curve, the stability scenario of the solutions is 
completely different in each case, as we demonstrate be¬ 
low. First we choose V = /R5- In both cases, a = ±1, 
no linearly unstable solutions were obtained, at least for 
moderately high amplitude solutions with P ~ 10 as is 
displayed in Fig. 3 and Fig. 4. 

With respect to the V = —/R5 case, we have found 
only unstable solutions in the a = 1 case as depicted 
in Fig. 5. The instability here appears due to an asym¬ 
metric internal mode that bifurcates from the continuous 
spectrum of (11), at A = ±b. This internal mode appears 
as soon as P > 0, thus there are no stable solutions be¬ 
cause this mode moves to the ImA > 0 region. Note that 
this instability rate is quite low when \b — bi\ <C 1. The 
Fig. 6 clearly shows that branches with higher values of 
7 tend to present a lower degree of instability in this case. 
In the V = — a/ 045 case we have found stable solutions 
only in the case of the defocusing nonlinearity, cr = — 1 
as is shown in Fig. 7. 


B. Effects of asymetry parameter /3 

The effect of different values of the asymmetry con¬ 
stant /3 has also been investigated. We have studied all 
four combinations of V = Tv/045 and a = ±1. See in 
Fig. 8 and Fig. 9 the case V = \/0d3- The effect of /3 in 
both a = ±1 cases is to shift b\ closer to b\ = 0.5 as ex¬ 
pected, with more discrete modes appearing as /3 — 1 
(See Fig. 2). The vertical asymptote b = b\ + <7/7 



v - 705, a = 0,7 = 10 , (7 = -1 



FIG. 9: Total energy P as a function of b with V — a/575, 
7 = 10 , <j — —1 and several values of saturation constant 7 . 
All solutions are stable. 
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FIG. 10: Total energy P as a function of b with V — — v / T5, 
7 = 10 , a = — 1 and various values of asymmetry constant ft. 
All solutions are stable. 


V = — v4l5, a = 0,7 = 10,a = 1 



FIG. 8 : Total energy P as a function of b with V — a/075, 
P = — 0.5, < 7=1 and different values of saturation constant 
7 . All solutions are stable. 

is still present in the branches. Regarding stability, no 
changes were observed: all solutions were found to be 
stable in the V = /R5 case as illustrated in Fig. 8 and 
Fig. 9. Solutions with V = — \/R5 are stable for <7 = — 1 
(See Fig. 10) and all solutions are unstable in the a = 1 
case (See Fig. 11). 

Note that in both cases, cr = =bl, with V = — 
there are stable solutions only if ft = 0, i.e., the TT- 
symmetric case. In Fig. 12 one may see that, in fact 
max{ImA} 0 as (3 0. 


FIG. 11: Total energy P as a function of b with V = — x/Odi, 
7 = 10, cr = 1 and different values of asymmetry constant /3. 
All solutions are stable. 


C. Investigation of gain and loss parameter a 

Here we fix f3 = —0.5, 7 = 10 and change a to inves¬ 
tigate its effects on the existence and stability of fun¬ 
damental branches. The first case to be investigated 
is the one with V = v^Odi and a = 1 and we turn to 
Fig. 13, from where one may conclude that all observed 
solutions are stable. The second case with V = \/0d3 
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P 


FIG. 12: Maximum value of the imaginary part of (11) as a 
function of total energy of localized modes P. Note that all 
(3 < 0 branches are unstable but max{ImA} — )► 0 as ft —>• 0. 
Parameters are V — 7 = 10, a = 1. 

and a = — 1 shows that an unstable region may exist 
above some threshold in P if a < 0 is close enough to 
a cr . The stable region becomes smaller as a gets closer 
to <a cr , eventually disappearing along with the branch it¬ 
self when a —)> a cr . Next, for V = — \/T5 and a = 1. 
we note that there are no stable solutions. Finally, for 
V = — \/T5 and a = —1, all solutions are stable. 



FIG. 13: Total energy P as a function of b with V = v / 0J5, 
7 = 10, a = 1 and various values of a. All solutions are 
stable. 



FIG. 14: Total energy P as a function of b with V — v / T5, 
7 = 10, a = — 1 and several values of a. Thick lines (lines) 
represent linearly stable (unstable) solutions. 



FIG. 15: Total energy P as a function of b with V = — v/075, 
7 = 10, a = 1 and different values of a. All solutions are 
unstable. 



0 0.05 0.1 0.15 0.2 0.25 0.3 

b 

FIG. 16: Total energy P as a function of b with V = — v / T5, 
7 = 10 , cr = 1 and different values of asymmetry constant ft. 
All solutions are stable. 


V. DYNAMICS OF LOCALIZED SOLUTIONS 

To corroborate the results described in the last section, 
in this section we proceed to make comparisons between 
the model of perturbations (11) and direct propagation 
of (3) using the solutions found as initial conditions in 
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(1). In all propagations we have added a 1% of noise to 
the amplitude at ( = 0. Here our objective is to show 
that the linear stability analysis predicts quite well the 
evolution of localized modes under weak perturbations. 


4 
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0 

-20 



FIG. 17: (Color online) (Left panel) Intensity \u\ 2 during prop¬ 
agation of a high power P — 10 stable soliton with b = 0.23. 
(Right panel) Corresponding eigenvalues of (11). Parameters 
are V — v / R5, 0 — —0.5, a = 0,7 = 10 and a — 1 . 


The first property that we have observed is the stability 
of three, among four possible combinations of signals of 
V and cr, three cases: V = y/&5 and a = =bl;R = — \/R5 
and g = — 1, exhibit stable solutions as possible to see in 
Fig. 17 and Fig. 18 for the V = \Z035 and a = 1 case. In 
all figures of dynamics we have used 7 = 10 as it does not 
change the existence or absence of stable solutions. Note 
that the discrete modes of stable solitons with /3 ^ 0 
tends to the negative imaginary part of the plane, dif¬ 
ferently from the (3 = 0 case, where the discrete modes 
of stable solitons must have ImA = 0 as it is depicted in 
Fig. 18 . 

In Fig. 19 we show two examples of unstable solutions 
in the branch with a = —0.14 (3 = —0.5. Note that 
while (3 does not change by itself the stability proper¬ 
ties of the solutions when all other parameters are fixed, 
the parameter a < 0 may in fact introduce instability in 
the system, both of oscillatory (complex A) and purely 
exponentially decaying (purely complex A) types, as dis¬ 
played in Fig. 19. Dynamics of unstable propagation in 
case V = — v/fhh and cr = 1 is shown Fig. 20. Finally in 
Fig. 21 illustrates the dynamics of a stable solution in 
the case V = —a/T5 and cr = —1. 

Looking at the stability spectrum of (11), one may con¬ 
clude that the role of the sign of V has with respect to 
the stability of branches may be understood by consid¬ 
ering that non-TT-symmetric cases, i.e., (3 ^ 0, do not 
possess the symmetry A = A* as expected. So, one may 
find discrete eigenvalues A in either ImA > 0 or ImA < 0 
regions. Now, if one inspects a branch where the po¬ 
tential is substituted as G G*, the spectrum of (11) 
changes to A —>■ A*, as all terms other than L are real in 
L s . Thus, whenever a stable solution exists with some 
ImA < 0 for a given G, there is an unstable solution in 


V = v^5, a = 0,7 = 10 , a = 1 


0=0 b = 0.1, P = 1.0498 



20 0 * -1 0 1 


Im A 

0 = -0.5 b = 0.1985, P = 1.0452 

2 ■ 



20 0 r -° 01 0 0 01 

4 ^ Im A 


FIG. 18: (Color online) (Left panels) Intensities \u \ 2 dur¬ 
ing propagation of stable solutions of branches with dif¬ 
ferent values of the asymetry constant /3. Parameters are 
V — v / T5, ol — 0,7 = 10 and a = 1 . 


the branch of the system with the potential G*. There 
are cases where changing G for G* does not change the 
stability of a function. This happens solely when a stable 
solution has ImA = 0 as possible to see in Fig. 17 and 
Fig. 21 for example. 
















6 = 0.0006367 



Im A 



t 20 0 C 


Im A 


FIG. 21: (Color online)(Left panel) Intensity \u\ 2 during prop¬ 
agation for a high power, P — 10, stable solution with b = 
0.0412. Parameters are V — —a/OA,/3 — — 0.5, = 0,7 = 10 

and a — — 1 . 


FIG. 19: (Color online)(Left panels) A Intensities \u \ 2 dur¬ 
ing propagation of unstable solutions. (Right panels) Cor¬ 
responding eigenvalues of (11) Note the that de linear sta¬ 
bility model correctly predicts the onset of a oscillatory in¬ 
stability in the solution with b m 0.0153 and the expo¬ 
nential decay in the case b = 0.0006367. Parameters are 
V = y/oH, P = -0.5, a = -0.14 ,7 = 10 and a = -1. 



FIG. 20: (Color online)(Left panel) Intensitt \u\ 2 during prop¬ 
agation of an unstable solution. (Right panel) Correspond¬ 
ing eigenvalues of (11). Note that the linear stability model 
correctly predicts the onset of a oscillatory instability. Also 
note that the instability develops quite slowly. Parameters 
are V — — v / R5, P — —0.5 ,a = 0,7 = 10 and am 1. 
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VI. CONCLUSIONS 

Summing up, we have studied the existence and stabil¬ 
ity properties of soliton families in media with a saturable 
nonlinearity in a non-TT-symmetric complex refractive 
index. We have found the unexpected result that such 
media does support the existence of asymmetric solitons 
in both self-focusing as well self-defocusing case. Fur¬ 
thermore, all branches of solitons found to bifurcate from 
the linear limit in case the saturation parameter 7^0, 
exhibit a vertical asymptote in the power-eigenvalue dia¬ 
gram, revealing that 7 restricts the range where localized 
modes are supported. Also, the asymptotic value of b for 
large powers is uniquely determined by the saturation pa¬ 
rameter and this value decreases for higher 7 driving the 
system to a less unstable situation. This result is similar 


to the result reported in TT-symmetric optical lattices [? 
]. Furthermore, as the amplitude of the soliton increases 
its profile becomes similar to the shape of the linear eigen¬ 
mode from where it was originated. In case of positive 

V we find stable solitons in both media, i.e., under the 
effect of self-focusing or self-defocusing. In cases where 

V takes on negative values we find that, when it comes 
to the stability properties these are complete altered and 
all solutions become unstable in the self-focusing media, 
i.e., for a = 1. 
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